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ABSTRACT 

We utilize high-resolution, high signal-to-noise spectra to perform a 
differential analysis of Fe abundances in the common proper-motion pair f 6 Cyg 
A and B. We confirm that both stars are slightly metal-rich compared to the 
Sun, and we show for the first time that the primary is enhanced in Fe relative 
to the secondary by a significant amount. We find A[Fe/H]= +0.025 ± 0.009. 
This tends to support the "self-pollution" scenario proposed by Gonzalez 
(1998), though lack of a complete understanding of small primordial metallicity 
variations among binaries and open cluster members prevents a definitive 
conclusion. 

Subject headings: stars: planetary systems — stars: individual (HD 186408, HD 
186427) - stars: abundances — stars:chemically peculiar 
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1. Introduction 

The nearby common proper-motion pair 16 Cyg A and B (HD 186408 and HD 186427, 
respectively) have been of particular interest to stellar astronomers for a number of reasons. 
In their physical characteristics, they are both very similar to the Sun and as such have been 
labeled "solar twins" or "solar analogs" (Friel et al. 1993, hereafter F93). Like any human 
twins, however, these two stars clearly cannot be exactly similar to the Sun or even to each 
other, and it is their differences that have generated the most recent interest. Most notably, 
radial-velocity studies have revealed that 16 Cyg B harbors a planetary-mass companion, 
while 16 Cyg A apparently does not (Cochran et al. 1997). This remarkable observation 
has itself been invoked by some to explain other unexpected chemical differences between 
these two stars (Cochran et al. 1997; King 1997; Gonzalez 1998, hereafter G98) - their 
Li abundances differ by a factor of five of more (Deliyannis et al. 2000, hereafter D00), 
while studies of the Fe abundance of this pair have consistently suggested that 16 Cyg A 
may be slightly more metal-rich (~0.05 dex) than its companion (G98). How are these 
differences achieved in such otherwise similar stars, which share a presumably common 
primordial environment? To what level can dynamic interactions with a planetary system 
affect abundances in a stellar photosphere? 

Unfortunately, with the exception of lithium, the small differences in line depths 
between these two stars combined with the relatively large uncertainties in abundance 
determinations to date have prevented any conclusive statements to be made about 
deviations between 16 Cyg A and B with regard to any other element, including iron 
(G98; D00); within the uncertainties, the measured chemical differences are all effectively 
zero. The primary goal of this work is to utilize a new differential line abundance analysis 
method to better constrain the relative [Fe/H] values between 16 Cyg A and B. Such a 
method, as shown by Langer et al. (1998), eliminates non-trivial uncertainties in oscillator 
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strengths, gf values. Further, we utilize a more extensive Fe linelist covering a wider range 
in lower excitation potential (xi) than previous published studies. In §2, we describe our 
observations and the data reduction, in §3 we describe our analysis methods and present 
results, and we discuss their significance in §4. 

2. Observations 

Our intentions from the start have been to perform a very careful spectroscopic analysis 
of 16 Cyg A and B, optimized for differential analysis. We determined at the outset to make 
every effort to minimize any systematic differences which might arise due to variations 
in our observations, reductions, and subsequent analyses of each star. With this goal in 
mind, we obtained optical spectra of 16 Cyg A and B within 30 minutes of each other on 
December 22, 1999, with the McDonald Observatory 2.7m telescope. We employed the 
coude echelle spectrograph (Tull et al. 1995) and a 2048x2048 Tektronix CCD, and made 
no changes to the instrument between the two observations. We achieved a resolving power 
of 58,000 (as measured on a Th-Ar lamp spectrum) with a signal-to-noise ratio (S/N) of 
~400 per pixel at 6700 A. Spectra of a hot star were also obtained within one hour and at a 
similar airmass in order to compensate for telluric features in the spectra of 16 Cyg A and 
B. 

Data for 16 Cyg A and B were reduced in as nearly an identical manner as possible, 
in both cases following the general method described in Gonzalez (1997). Each image 
was processed with the same bias and flat fields, and corrections for scattered light were 
made with the same fitting functions in both spectra. In addition, continua of individual 
orders were normalized with exactly the same functions in both stars. These steps were 
undertaken to minimize possible systematic differences between line depths which would 
ultimately impact our differential analysis. Our efforts yielded continuous one-dimensional 
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spectra in the blue up to 5500 A, with gaps thereafter up to roughly 10,000 A. 

3. Analysis 

3.1. Equivalent Widths 

To eliminate potential bias in measuring the equivalent widths (EWs) of spectral 
features, one of us (GG) renamed each spectrum, and the other (CL) performed the 
measurements and analysis with no prior knowledge of which spectrum was associated with 
which star. We employed a list of high-quality atomic lines used in our previous studies of 
planet bearing stars, and measured the EWs of 60 Fe I and 8 Fe II lines spanning broad 
ranges in xi from 0.09 to 6.22 eV and in EW from 4.7 to 119.8 mA. In every case the EW 
was measured for a line first in one star, then immediately afterwards in the second star 
before moving on to the next line. This procedure was followed to minimize any systematic 
differences in EW measurements between the two stars. Telluric contamination in some 
Fe lines was addressed by dividing the science object spectra by that of a hot star, and 
where performed, this division, as well as considerations of blending by nearby features, 
was treated identically for each line in each star. Our linelist and measured EW values are 
presented in Table 1. 

3.2. Standard Analysis 

We employed the measured EWs from §3.1 and associated gf values taken from 
Gonzalez et al. (2000) to derive a set of basic stellar parameters (T eS , \ogg, £ t , and 
[Fe/H]) and associated uncertainties for 16 Cyg A and B, using the methods described 
extensively in G98 and Gonzalez & Vanture (1998). Briefly, all four of these parameters 
were systematically iterated with a recent version of the LTE abundance code MOOG 
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(Sneden, 1973) using Kurucz model atmospheres (1993) until the mean Fe I and Fe II 
abundances were equal and correlations of the individual Fe I line abundances with both xi 
and the logarithm of reduced equivalent width (REW) were zero. 

Following the method described in Gonzalez & Vanture (1998), uncertainties in £ t 
were determined from the standard deviation in the slope of a least-squares fit to the Fe I 
vs. \og(REW) data. This resulting uncertainty in £ t was summed in quadrature with the 
standard deviation in the slope of a least-squares fit to the Fe I vs. xi data to estimate 
the uncertainty in T eS . We combined theses uncertainties in x\ an d T e g with the observed 
scatter in modeled Fe I line abundances to calculate our final uncertainty in [Fe/H].[| 

Our results using these standard methods are presented in Table 2, and we confirm 
previous studies: 16 Cyg A and B are very similar to each other, as well as to the Sun. The 
pair are each slightly metal-enhanced relative to solar, however, with 16 Cyg A presenting 
an 0.03 dex overabundance of Fe relative to 16 Cyg B, although this difference is not 
significant given our formal estimates of uncertainty. 

3.3. Differential Analysis 

To reduce the uncertainties in our estimates from §3.2, we reanalyzed our EW data 
from these two stars employing a technique similar to that described and effectively 
utilized by Langer et al. (1998) in searching for small metallicity variations amongst red 
giants in M92. Essentially, this method mirrors the "standard" method discussed above, 
but instead of determining chemical abundances for each star by averaging the results 
determined from individual lines and subsequently comparing these averages between stars, 
it determines differential chemical abundances between stars by taking the differences 

throughout this paper, all stated uncertainties are at the 1-cr level. 
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in abundances calculated for each line individually. This differential strategy effectively 
eliminates uncertainties in gf values, and allows one to take advantage of the increased 
precision in A[Fe/H] to further constrain differential values of the other stellar parameters. 

We used our values for T eS , log g, and £ t from §3.2 to calculate [Fe/H] for each 
line in each star, then determined differential Fe I abundances (A [Fe/H]) for each line 
individually.0 If the AT e fj and A£ t values from §3.2 were both correct, then plots of A [Fe/H] 
vs. both x\ an d log(REW) would be uncorrelated. Our initial such plots using these "A" 
values as inputs, however, showed slopes of 0.003 ± 0.002 dex/eV for A [Fe/H] vs. xi an d 
-0.010 ±0.009 dex/\og(REW) for A [Fe/H] vs. log(REW) (for comparison, these values in 
our standard analysis of 16 Cyg A in §3.2 were 0.001±0.006 and 0.003±0.024, respectively). 
While these results are only slightly larger than their standard uncertainties, we note that 
they are both more precise and appreciably larger in magnitude than the final "zero-slope" 
conditions which resulted from our standard analysis. Such correlations indicate that the 
values of A T e g and A£ t calculated from the results in Table 2 are off by some amount, and 
so we iterated AT fj and A£ t by changing T e ff(i 6 c y gB) and £t(i6CygB) until these correlations 
were forced to zero.0 Figure 1 presents a visualization of the effect of changing AT e g and 
A£t on the x\ an d log(REW) slopes. For logg = 4.21 and 4.26 for 16 Cyg A and B, 
respectively, we find our best solution at AT cff = 62 ± 14K and A£ t = 0.05 ± 0.01 km/s, 
where the stated uncertainties were calculated in precisely the same manner as described 
in §3.2. We note that these values are in very good agreement with our earlier, though less 
precise, results generated via the "standard" method. 

Figures 2 (a) and (b) show our final plots of A [Fe/H] vs. both xi an d log(REW) for 
our optimal values of AT e g and A£ t ; immediately apparent is the fact that for the vast 

2 In the following discussion, all "A" values are of the form AX = A 16Cyg A — A 16 c yg B- 
3 Additionally, one consistently discrepant Fe I line at 6864.32 A was discarded. 
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majority of lines analyzed, A[Fe/H] is positive. We find a mean value of A[Fe/H] and its 
associated uncertainty of +0.025 ± 0.009 dex, where the stated uncertainty in A[Fe/H] was 
calculated via the method described in §3.2, utilizing the uncertainties in AT c fj, A£ t , and 
the observed scatter of ±0.021 dex in individual line A[Fe/H] values. Uncertainty in Alogg 
produced an effect of ±0.001 dex in A[Fe/H], and was not considered in the stated error of 
the mean value of A[Fe/H]. This uncertainty is in fact dominated by uncertainty in AT e fr, 
with a ±14K change effecting a ±0.008 dex change in A[Fe/H]. 

4. Discussion 
4.1. Reliability 

We note that the atmospheric parameters for 16 Cyg A and B calculated using the 
"standard" method stand in good agreement with previous published studies - a heartening 
fact in light of discrepancies in EW measurements^ as well as variations in methodology 
among the various authors. We believe our absolute values are more reliable than those of 
other recent studies due to our use of more Fe I lines spanning a larger range in \\ an d EW 
(D00, for example, only employed 20 Fe I lines with the smallest xi an d EW values being 
2.18 eV and 24 mA, respectively). Furthermore, using our method of differential analysis, 
we are able to fine-tune the observed values of AT e g- and A£ t , and as a result have shown 
for the first time a statistically significant difference in [Fe/H] between 16 Cyg A and B. In 
Figure 3 we compare our resulting A [Fe/H] values with those of previous studies. 

4 Two typos in G98 were discovered in the course of preparing this paper: EWs for the Fe 
I feature at 6710.32 A were mistakenly repeated for 6733.15 A (the latter was not measured 
in that study), and EWs for the Ti I feature noted at 6126.22 A were actually from a Ti I 
line at 6261.11 A. 
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The magnitude of this difference remains quite small, though, and while we have 
strived to eliminate possible systematic differences in our observations and analyses of these 
two stars, we must consider other potential mechanisms which might mimic true differences 
in chemical abundance before we can conclude that 16 Cyg A is in fact enhanced in Fe 
relative to 16 Cyg B. Most notably, measured values of [Fe/H] are known to vary with 
chromospheric activity levels. A recent study of Ca II fluxes in planet-bearing stars by 
Henry et al. (2000), however, shows low values of i?^K f° r 16 Cyg B, indicating little 
variability and surface activity. Our own spectra confirm the low levels of Ca II H and 
K emission, and show no signs of significant differences between 16 Cyg A and B in this 
regard. Furthermore, Hipparcos data sets the photometric variability of each of these stars 
at the 0.0007 magnitude level. We therefore are confident that differences in chromospheric 
activity do not play a significant role in accounting for the measured values of A [Fe/H]. 

In addition, we examined possible systematic errors which might be introduced due to 
assumptions in our atmospheric modeling, such as the possibility that a difference in the 
temperature minimum between the two stars might yield differing values of [Fe/H] for the 
same T cff (Wallerstein, 1972). To test against this, we compared the central absorption in 
the Mg I triplet near 5170 A, and found no evidence supporting differences in temperature 
minima between 16 Cyg A and B. We further argue that given the common ages of 9 ± 2 
Gyr (G98) and similar values of T eS , logg, and luminosity of these two stars - confirmed 
independently from photometric analyses (c.f. F93) - it is unlikely that any systematic 
errors in this strictly differential study might be due to more subtle model assumptions, 
such as that of LTE. We therefore conclude that our measured value of A [Fe/H] is reliable 
to the level of our stated uncertainty, and that the photosphere of 16 Cyg A is significantly 
more iron-rich than that of 16 Cyg B. 



4.2. Possible Explanations 



An obvious explanation for the observed value of A[Fe/H] is that it represents a 
primordial difference in the chemical composition of 16 Cyg A and B. Little high- resolution 
work exists on abundance differences between members of multiple star systems, however, 
and models of the formation of stellar systems are currently not sophisticated enough to 
address deviations in chemical abundance at the level of precision that we report here. We 
intend to actively pursue the former shortcoming by employing our differential method on 
a wide sample of stellar systems; however, until such primordial fluctuations are at least 
empirically, if not theoretically constrained, there is little more that can be said for or 
against this hypothesis. 

If, on the other hand, the observed difference in [Fe/H] is not entirely primordial, then 
it must be set by variations in the evolutionary history of these two stars. A potential clue 
to the responsible mechanisms may lie in the aforementioned fact that the Li abundances of 
16 Cyg A and B differ by a factor of at least five (D00). While Ryan (2000) has noted that 
planet-bearing stars do not as a whole appear to differ from the field population in their 
values of [Li/H] - a finding supported by the larger sample discussed in Gonzalez et al. 
(2000) - it remains a challenge for standard theories of stellar atmospheres and evolution 
to explain the large difference in [Li/H] between 16 Cyg A and B given the considerable 
similarity of their physical characteristics and presumably shared environmental history. 

In response, various authors have discussed the possibility that the chemical composition 
of stellar photospheres can be affected by the presence of planetary companions. Cochran 
et al. (1997) and King (1997) proposed that dynamical interactions between a rotating star 
and its protoplanetary disk might significantly alter the angular momentum evolution of 
the star, and hence the rate at which the star depletes Li. This model is well supported by 
the recent investigations of Li and Be abundances reported by D00, wherein these authors 
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conclude that models of slow rotational mixing can explain the relative abundances of Li 
and Be in 16 Cyg A and other anomolously Li-rich F and G stars. Such a mechanism 
cannot explain the more subtle variation in Fe which we report here, however. 

Alternately, G98 proposed that the increased lithium content of 16 Cyg A's photosphere 
may be the result of that star having consumed planetary material in its outer convection 
zone. Such self-pollution by materials of either chondritic or gas-giant composition would 
indeed produce an increase in the Li abundance, while only slightly affecting the Fe and 
Be abundances (G98, D00). If our result indicating that 16 Cyg A is slightly enhanced in 
Fe relative to 16 Cyg B is not a primordial effect, then it lends tentative support to this 
scenario. From the estimates presented by G98, we calculate that the observed difference in 
[Fe/H] between 16 Cyg A and B can be explained by the accretion of 2.5 M ffi of chondritic 
or 0.3 Mj of gas giant material. The amount of Li enhancement due to accretion of these 
bodies is difficult to determine, given the marked non-linear time, T e e, and dynamical 
dependence of Li depletion, along with the unkown timing of the putative accretion event (s). 
Nevertheless, it remains plausible that the observed difference in Li between 16 Cyg A and 
B is explainable by accretion at some intermediate age. 

The authors are grateful to George Wallerstein for helpful discussions and support 
through the Kenilworth Fund of the New York Community Trust, and to the anonymous 
referee whose helpful comments greatly increased the clarity of our manuscript. This work 
has utilized the Simbad database at CDS, Strasbourg, France, in addition to the abstract 
services of ADS. 
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Table 1. Fe I and II Equivalent Widths - 16 Cyg A and B 
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Table 1 — Continued 
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Table 2. Spectroscopically Determined Physical Parameters of 16 Cyg A and B using 

"standard" method. 



Parameter 16 Cyg A 16 Cyg B Sun 
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Fig. 1.— A plot of the slopes of log(REW) (=log(EW/A)) and xi vs. A[Fe/H] for a range 
of AT eff and A£ t . For each series, \ogg = 4.21 and 4.26 for 16 Cyg A and B, respectively. 
Each point represents a change of 5 K in temperature difference, from AT c ff = 80 K for the 
leftmost points in each series. Shown with a cross is the (0,0) coordinate representing the 
formal solution; our optimal solution is AT eff = 62 K and A£ t = 0.05 km/s. 

Fig. 2.— A[Fe/H] vs. (a) x\ and (b) log(REW) for our optimal values of AT cff , Af t and 
A log g. The solid lines are least-squares fits to the data. 

Fig. 3. — A comparison of A[Fe/H] values for 16 Cyg A and B from various published studies, 
with one-sigma error bars. 
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